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Resumen 


En este documento se presenta un procedimiento para determinar las 
constantes empíricas de von Kármán « y de adición A de la ley de 
distribución de velocidad logarítmica, y se fundamenta en una 
estimación global de los valores mínimos de error entre el modelo de la 
ley logarítmica y los datos de una medición. Para verificar la bondad del 
método propuesto, se realizó una comparación con otros métodos 
presentados en la literatura, utilizando los datos de un perfil de 
velocidad obtenidos en un canal de pendiente variable de sección 
rectangular, con un velocímetro acústico de efecto Doppler (ADV, 
Acoustic Doppler Velocimeter). La prueba experimental fue en flujo 
uniforme con fondo de rugosidad ks = 0.000535 m y un número de 
Reynolds R. > 36 000. El método propuesto permite estimar de forma 
simultánea las constantes k y A, localizando un valor óptimo 
denominado Emi en un plano bidimensional de estimación de error de 
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las constantes. Los resultados del método propuesto en comparación 
con otros reportados en la literatura indican una mejora o al menos una 
estimación similar. Además, se comprobó el ajuste de la ley logarítmica 
con los valores obtenidos de la prueba experimental en el rango 
establecido 30 < y+ < 0.38+, utilizando los valores calculados de « y A con 
el método propuesto. 


Palabras clave: ley logarítmica, constante de von Kármán, constante 
de adición, velocimetría acústica Doppler. 


Abstract 


This paper presents a procedure to determine the empirical constants of 
von Kármán kx and additive A of the law of logarithmic velocity 
distribution, is based on a global estimation of the minimum values of 
error between the logarithmic law model and the data of a 
measurement. To verify the goodness of the proposed method, a 
comparison was made with other methods presented in the literature, 
using the data of a velocity profile obtained in a variable slope channel 
of rectangular section, with an Acoustic Doppler Velocimeter (ADV). The 
experimental test was in uniform flow with bed roughness ks= 0.000535 
m and Reynolds number Re. > 36 000. The proposed method allows to 
simultaneously estimate the empirical constants x and A, locating an 
optimal value called Emin in a two-dimensional plane of error estimation 
of the constants. The results of the proposed method compared to 
others reported in the literature indicate an improvement or at least a 
similar estimate. In addition, the adjustment of the logarithmic law with 
the values obtained from the experimental test in the established range 
30 < y+ < 0.39+ was verified, using the calculated values of « and A with 
the proposed method. 


Keywords: Logarithmic law, von Kármán constant, additive constant, 
acoustic Doppler velocimetry. 
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Introducción 


La distribución de velocidades medias en un cauce natural o artificial se 
considera por simplicidad como un proceso bidimensional que tiene 
regiones con diferentes características de velocidad en función del 
tirante (Nezu 8 Nakagawa, 1993). Estas regiones se evalúan por la 
relación y/h, donde y es la coordenada vertical y h es la profundidad o 
tirante del flujo. En la región cercana al fondo 0 < y/h < 0.2, la velocidad 
principal es descrita por la clásica “ley de pared” (Coles, 1956), que se 
compone de tres zonas: subcapa viscosa, amortiguamiento y logarítmica 
(Petrie € Diplas, 2016; Spalding, 1961). La existencia de la subcapa 
viscosa y la zona de amortiguamiento está en función de la rugosidad 
del fondo, por ejemplo, si el valor de escala de rugosidad adimensional 
k? <5 (k? = ksu,/v, donde k, es la rugosidad absoluta, u, es la velocidad 
de corte y v es la viscosidad cinemática), el flujo se considera con fondo 
hidráulicamente liso y están presentes las tres zonas. Si 5<k¿ <70, el 
flujo penetra la zona de amortiguamiento y la subcapa viscosa 
desaparece, considerándose como flujo con fondo incompletamente 
rugoso. Finalmente, si ki >70, se considera como flujo con fondo 
completamente rugoso y las dos primeras zonas desaparecen, 
permaneciendo sólo la zona logarítmica. En el rango 0.2 < y/h<1.0 el 
perfil de velocidad no sigue una tendencia logarítmica, ya que se tienen 
dos comportamientos: el denominado fenómeno de ¡inmersión 
(Bonakdari, Larrarte, Lassabatere, € Joannis, 2008; Yang, 2005) o de 
estela (Yang, 2009), donde los valores de una medición difieren 
respecto a los estimados con la “ley logarítmica”, ecuación (1). 


ut =Iny*+A4 (1) 


Donde u* =u/u, es la escala de velocidad adimensional; y* = yu,/v, la 
escala de longitud adimensional; u, la velocidad media en dirección 
principal del flujo; x, la constante de von Kármán, y 4 es la constante de 
adición. Este modelo logarítmico se ha utilizado para caracterizar el 
campo de velocidad de flujos turbulentos en ríos (Ferro 8 Baiamonte, 


1994; Petrie € Diplas, 2016), canales (Auel, Albayrak, € Boes, 2014; 
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Tominaga € Nezu, 1992) y tuberías (Bailey, Vallikivi, Hultmark, 8 Smits, 
2014; Perry, Hafez, 8 Chong, 2001). 


De acuerdo con Segalini, Órlú y Alfredsson (2013), en los trabajos 
experimentales que se exponen en la literatura, el rango donde la ley 
logarítmica reproduce de modo correcto el perfil de velocidad no está 
completamente definido. Según Segalini et al. (2013), el límite inferior 
oscila entre 30< y* < 300 e incluso hay publicaciones que indican un 
límite hasta 600 (McKeon, Li, Jiang, Morrison, € Smits, 2004; Zagarola 
8 Smits, 1998); para el límite superior se encuentra alrededor de 
0.1<y*/8* <0.3 (6* =8u,/v, donde Y denota el espesor de la capa 
límite); también hay publicaciones que muestran valores hasta 0.75 
(Zanoun, Durst, € Nagib, 2003). 


Para el uso de la ley logarítmica es necesario estimar los valores de las 
constantes empíricas k y A, y existen diferentes métodos para 
determinarlas en diversas condiciones de flujo, a partir de datos de una 
medición (Monkewitz, 2017; Segalini et a/., 2013; Zanoun et al., 2003). 
Los parámetros básicos que deben tomarse en cuenta para estimar las 
constantes son número de Reynolds y Froude (Bailey et al., 2014; 
Monkewitz, 2017); rugosidad del fondo (Auel et al., 2014), y la relación 
de aspecto del cauce b/h (en el caso de una sección rectangular, donde 
b es el ancho del fondo del canal) (Bonakdari et a/., 2008). 


En este documento se presenta un procedimiento para determinar las 
constantes empíricas xk y A de la ley logarítmica, y se fundamenta en 
una estimación global de los valores mínimos de error entre el modelo 
de la ley logarítmica y los datos de una medición. Además, con el fin de 
conocer los alcances de este método, se realizó una comparación con 
otros métodos descritos por Segalini et al. (2013), que son mínimos 
cuadrados, función de indicador y dispersión xk-—A. En el apartado 
siguiente se describen tales métodos, junto con el propuesto en este 
trabajo, denominado topología. 


La determinación de las constantes empíricas ante diferentes valores de 
número de Reynolds y Froude, rugosidad del fondo y relación de aspecto 
permite conocer sus variaciones y, por tanto, es posible utilizar con 
mayor precisión metodologías para determinar la velocidad de corte 
mediante el uso de la ley logarítmica, como la expuesta por Mendoza- 
González y Aguilar-Chávez (2018), donde se asumen valores fijos para xk 
y A. 
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Métodos de estimación de las constantes rx y A 


Mínimos cuadrados 


En el método de mínimos cuadrados se estima la constante de von 
Kármán con la ecuación (2) (Segalini et a/., 2013). 


N Buln? y? (2imy ] 
a» + [L,...,N) (2) 


Donde N es el número de valores de una medición que se consideran 
dentro del rango de la ley logarítmica e ¿ denota el contador de cada 
valor. Una vez calculado x, se procede con la estimación de A mediante 
la ecuación (3). 


A=2 Yi (ut —Inyf) (3) 


Función de indicador 


Este método tiene su origen en aplicar la derivada a la escala de 
velocidad u* respecto a la escala de longitud y* de la ley logarítmica. 
Entonces, se tiene la siguiente relación, ecuación (4). 
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2 = (Eny* +A) (4) 
dy+  dy* Ax y 


Desarrollando la derivada se tiene la ecuación (5). 


du* 


11 
> (5) 


Ky 


Finalmente, despejando la inversa de x, se obtiene la función de 
indicador, ecuación (6). 


Le yo (6) 


rg dy+ 


Para aplicar la ecuación anterior a valores discretos de una medición se 
utiliza la ecuación (7). 


INE! + du* 
E 
Ko oN ¿Yi dy+ 


(7) 


i 


Mediante la ecuación (7) se calcula xk, mientras que A se determina con 
la ecuación (3). Este método también ha sido utilizado para discernir el 
rango de la región logarítmica. Este procedimiento es expuesto por 
Zanoun et al. (2003). La dificultad de tal método es determinar la 
derivada del perfil de velocidad du*/dy*, que tiene una gran dispersión 
debido a la incertidumbre inherente de los datos de una medición, 
siendo mayor en la cercanía de la pared (Segalini et a/., 2013). Por lo 
tanto, para aplicarlo método se recomienda tomar las muestras de 
velocidad a una distancia constante (Ósterlund, 1999) y emplear un 
esquema de diferenciación central de orden superior (Buschmann 8 
Gad-El-Hak, 2003). 


Dispersión k-— A 
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Este método fue propuesto por Alfredsson, Imayama, Lingwood, Orlú y 
Segalini (2013), y consiste en un reordenamiento de la ley logarítmica, 
tal como lo indica la ecuación (8). 


KA = ku* — Iny?* (8) 


Para emplear la ecuación (8) se propone un rango de M valores de xk, y 
este rango se divide en ¡€ [1,...,Mj valores de evaluación, para los N 
valores de la región logarítmica. Con ello se obtiene una distancia de 
error (AxA;), para cada valor de x propuesto, tal distancia se define por 
la diferencia entre el max(xA,) y el min(xA;,). El valor óptimo de ajuste 
para x, es aquel donde se tiene el min(AxA;);. Por último, A se calcula 
con la ecuación (3), o bien, a partir de los valores xA,, con la ecuación 


(9). 


5 — Y: KA; (9) 


Topología (propuesta) 


El método propuesto consiste en un reordenamiento de la ley 
logarítmica, como se muestra en la ecuación (10). 


u*—-Iny*-A=0 (10) 


En la expresión anterior, los valores de una medición u* y y* deben 
guiar el resultado a cero; en un caso ideal se podría tener un solo valor 
de x y A, para indicar el valor exacto del perfil de velocidades de una 
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medición. En este caso, se propone que existe una pequeña diferencia 
entre los valores de una medición y el valor de estimación de la 
ecuación (10), y se define como error de estimación E;. La forma 
funcional del error se aplica sobre un espacio paramétrico de los valores 
de x y A, por lo que E;: Q(x, A) = [Kin Kmax] X Amin Amar], donde M(x, A) es 
un espacio métrico para determinar los valores de las constantes del 
modelo logarítmico. Entonces, por cada valor de una medición se 
obtiene un valor de error sobre el espacio £M, tal como se indica en la 
ecuación (11). 


E¡(x, A) = u? - In y¿ —A (11) 


Con los valores de E;(x, 4) se obtiene el plano métrico del error, que en 
teoría se considera continuo; sin embargo, la forma funcional de x y A 
no son conocidas, por lo que se considera una evaluación discreta por 
valores conocidos de éstas. Para obtener una medida escalar de la 
magnitud del vector de error de estimación E, se aplica alguna norma de 
Hólder, ecuación (12). La selección de una norma se discute después. 


E = [E,(x, A) ll, = (ln ] pauta 00 (12) 


Donde ll-| es el operador norma y p € [1,0] indica el número de norma. 
De acuerdo con la definición formal de norma vectorial, la ecuación (12) 
puede ser cero, si y sólo si E;(x, 4) = 0 (Bernstein, 2009); esto indicaría 
un ajuste exacto de la ley logarítmica con los valores u* de una 
medición. En la práctica no es posible tal condición y, por tanto, ninguna 
combinación de valores kx y A permite igualar con cero la ecuación (12), 
sin embargo, existe una combinación que minimiza el error de 
estimación E y genera un ajuste adecuado de la ley logarítmica con un 
perfil de velocidad medido; entonces se dice que E;(x,4) -> 0 para una 
combinación de valores de x y A, en el dominio M(x, 4). 


En forma práctica, para estimar los valores de x y A que minimizan E, se 
propone una discretización del espacio del dominio, de modo que se 
tienen los siguientes vectores: x;=(1X,,..,Km) Y A¡=(41..,40), de 
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dimensión M y Q, respectivamente, además M,Q>0 y M,Q€Z. El 
espaciamiento Ak = (Kmax — Kmin)/M Y AA = (Amax — Amin)/Q pueden ser 
establecidos en función de la precisión que se requiera en la estimación 
de x y 4. Para establecer el rango de [Kinin max] Y [Amin Amax] se utilizan 
los valores reportados para estas constantes por otros autores, como los 
mostrados en la Figura 1. 


6.6 
5.8 
5.0 x= 
4.2 
3.4 


OHW0W0OMO00NANNDONATY0ODAMTIOOMNTRA 
LN LON LN LN WWWWRO0OIMIMAACAOANSMODOO 
DAMSDODODODODODODODODODODODODOODONOAANMDMOOOOOO0OO 
AA AAA AA AA AA AA AAA ANNNUNUNANAIUNN 
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5 DU=0 TEA SE 44 yw 
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zZ E E (0) INXGOVOVOY 2'm= 
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Figura 1. Valores de las constantes de la ley logarítmica obtenidas de 
varias investigaciones. 
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Una vez conocido el rango de combinación de los valores x; y A, se 
obtiene una matriz de error E(x;, Aj), donde el valor óptimo de x; y A, es 
aquel que min||£¿(«;, A,)!|,. 


Prueba experimental 


A fin de comparar los métodos de estimación de x y A, se llevó a cabo 
una prueba experimental en el canal de agua de pendiente variable del 
Instituto Mexicano de Tecnología del Agua (IMTA), que tiene geometría 
rectangular de ancho b = 0.245 m y una longitud L=5 m (ver Figura 2). 
La estación experimental cuenta con una bomba de 10 hp que abastece 
el caudal; un vertedor de aforo calibrado con la norma ISO 1438:2008, 
y una válvula para la regulación del caudal, entre otras componentes 
particulares, que permiten condiciones estables para los experimentos. 


O mt << 
! AM 


A "uu Pm 


AI AAA 


Figura 2. Estación experimental. 


Se implementó la técnica de velocimetría acústica de efecto doppler 
(ADV, Acoustic Doppler Velocimetry) para la medición de la velocidad 
del flujo, utilizando un equipo Vectrino Profiler”" de Nortek) (ver Figura 
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3). El equipo se posicionó a 3.5 m de distancia de la entrada del flujo 
para evitar defectos en el perfil de velocidad causados por la entrada del 
flujo o por la salida de éste en caída libre al tanque de recirculación. La 
frecuencia de muestreo fue de 100 Hz, con un tiempo de 
almacenamiento t=60 s. Como parte del experimento se instaló un 
tapete de lija en el fondo del canal, con rugosidad absoluta k, = 0.000535 
m (ver Figura 3). Las características del flujo y el valor de k, dieron 
como resultado la condición de flujo con fondo incompletamente rugoso 
(ver Tabla 1). 


Figura 3. Fondo de lija y equipo ADV Vectrino Profiler””, 


Tabla 1. Condiciones experimentales. 


Descripción Unidades Valor 
Pendiente S m/m 6.380 (10%) 
Ancho de base del canal b m 0.245 
Tirante h m 0.156 
Relación de aspecto b/h m/m 1.571 
Velocidad de corte u, m/s 2.062 (10?) 
Temperatura T eE 30.000 
Viscosidad cinemática v m?/s 8.007 (107?) 
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Caudal Q m*/s 1.631 (107?) 
Velocidad media U = Q/A m/s 0.427 
Número de Froude F. = U//gh (-) 0.345 
Número de Reynolds R, = UR, /v a 

e e (-) 36.593 (105%) 
R, = bh/(b + 2h), radio hidráulico 
Rugosidad adimensional k; (-) 13.777 


En el experimento se estableció la condición de flujo uniforme con una 
pendiente del canal S = 0.000638 (ver Figura 4). 


Figura 4. Condición de flujo uniforme en prueba experimental. 


La velocidad de corte se calculó con el modelo de Guo y Julien (2005), 
ecuación (13). 


=p 1) o Das] a) 


Donde gy es la aceleración de la gravedad; S, la pendiente del fondo 
(considerando flujo uniforme), y b es el ancho de base del canal. Este 
modelo contiene un factor de corrección que disminuye el error de 
estimar directamente la velocidad de corte con el modelo clásico 
u. = /ghS, sobre todo en bajas relaciones de aspecto b/h < 10. 


La posición y de los puntos de medición se obtuvo de forma directa con 
el equipo ADV, que permite medir la distancia al fondo y ubicar la 
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posición de las celdas de medición. Además, el ADV cuenta con un 
termistor para medir la temperatura del agua, a partir de la cual se 
asignó un valor de viscosidad cinemática, considerando las tablas de 
Wagner y Kretzschmar (2008), que están basadas en la formulación 
industrial ITAPWS-IF97. A manera de resumen, se presenta la Tabla 1 
con las condiciones experimentales de la prueba. 


Resultados 


Los métodos de estimación se aplicaron en el rango 30 < y* < 0.38* de 
los valores experimentales, donde d$* = 1558 (considerando un valor de 
espesor de capa límite $ = 6.05 cm, que se estimó en función del máximo 
valor de velocidad del perfil de velocidad). Este rango se plantea de 
acuerdo con los límites descritos por Segalini et a/. (2013), mismos que 
se mencionan en la introducción de este documento. 


En la metodología de análisis topológico se requieren establecer 
inicialmente los vectores x, y A,; en este caso se propone que xk € 
[0.30,0.45] y A € [3.40,6.60], de acuerdo con los rangos de la Figura 1. El 
espaciamiento Ax y AA de los valores de cada vector x; y A, está en 
función de la precisión que se quiera obtener, como ya se indicó. 


En la Figura 5 se muestra la topología de errores obtenida al evaluar la 
ecuación (12) con los valores de la prueba experimental, utilizando la 
norma p = 1 y un espaciamiento Ax = AA = 0.001; los valores de x y A que 
mejor ajustan el modelo teórico con los datos experimentales fueron 
0.418 y 5.067, respectivamente. 
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45.00 
18.98 
8.00 
3.37 
1.42 
0.450 0.60 


Figura 5. Topología de errores E, para p= 1 y Ax = AA = 0.001. 


En la Figura 6 se presentan los valores de min]|£;(x;, 4,)|], para diferentes 
valores de norma p y de aproximación Ax,A4. 
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meo p=1 -:e-p=2 -:“-p=5 --w-p=10 
--9-p=50 --e- p=100 --b-p=00 


0,1 0.01 1E=3 1E-4 1E=5 1-6 
Ax, AA 
Figura 6. Valores de Ejpin = min]|E, (x;, 44)!|, en función de Ax y AA para 
diferentes valores de p. 


De acuerdo con el resultado presentado en la Figura 6, a partir de un 
espaciamiento Ax =AA=0.001, el Ein se mantiene aproximadamente 
constante para cualquier valor de norma; esto indica que el aumento en 
la resolución de los vectores (Ax,AA4 < 0.001) no produce una disminución 
del error. 


Para la prueba experimental se determinaron los valores de x y A por los 
métodos de mínimos cuadrados, función de indicador, dispersión xk—A y 
topología. Los resultados se muestran en la Tabla 2. Para el método de 
topología se utilizó un espaciamiento Ak=AA=0.001; en el caso del 
método de dispersión xk — A, se utilizó la misma resolución para el vector 
de Kj. 


Tabla 2. Valores de x y A obtenidos de cada método. 


Método K A 
Mínimos cuadrados 0.425 5.308 
Función de indicador 0.385 4.114 
Dispersión k — A 0.431 5.469 
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p= 0.418 5.067 
p=2 0.425 5.308 
p=5 0.434 5.557 
Topología p=10 0.434 5.565 
p=50 0.431 5.503 
p=100 0.431 5.503 
p=00 0.431 5.504 


Con el fin de hacer una comparación entre los diferentes métodos se 
proponen los siguientes criterios de estimación del error: a) error 
absoluto acumulado (E,.,) ecuación (14), y b) error cuadrático 
acumulado (E..), ecuación (15), (Wackerly, Mendenhall, 8 Scheaffer, 
2007). 


+ + 
= Lu; . — UU; Se 
aa y | l experimental l teórico 


E 


(14) 


+ + 
= lu: , —Uu; dEi 
ca Zi | l experimental l teórico 


E (15) 


En la Figura 7 se presenta la comparación del E, Y Es. en la estimación 
de las constantes k y A por los diferentes métodos expuestos en este 
documento. 
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Ea a 


a a] X= HH (NN LON o o o o] 
a] LN 
hS he) Y 1 1 1 a 1 
o ln] E Il Il 
5 3 a, a a, Il Q 
pz 5 == a a 
a) E (7) a, 
3 3 
o , 
ts a E Topología 
2 c a 
E 59 
£ fa 
= 5 
E 


Figura 7. Valores de E,, Y Esa para cada método de estimación. 


De acuerdo con la Figura 7, el método con menor valor de E,, es el de 
topología con p = 1, mientras que el menor valor de E, se obtuvo con el 
método de mínimos cuadrados y con el de topología con p=2. El 
método de función de indicador contiene los más altos valores de error, 
esto es generado por el término diferencial de la ecuación (7), a pesar 
de que se aplicó una discretización centrada de segundo orden. En el 
caso del método de dispersión xk — A, sus valores de error no difieren de 
forma significativa de los valores obtenidos con el método de mínimos 
cuadrados, sin embargo, en ningún caso de forma de evaluación del 
error optimiza el ajuste. En el método de topología, si se consideran 
valores de p > 2 no se optimiza el ajuste en ninguno de los dos casos de 
evaluación del error. 
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El perfil adimensional de velocidades medias con los valores de las 
constantes xk y A obtenidas por el método de topología con la norma 
p=1 se presenta en la Figura 8. En este perfil se observa el fenómeno 
de estela e inmersión para y* > 400, donde los valores experimentales 
están por encima y debajo de la ley logarítmica estimada; en esta zona 
no se aplicaron los criterios de estimación de error. 


25 
20 
5 15 
10 E] Rango 30 < y* < 0.3 5* 
Ley logarítmica, 
k«=0.418 , A=5.067 
O Valores experimentales 
5 
101 10" y 10* 10* 


Figura 8. Perfil adimensional de velocidades medias. 


En este trabajo no se analiza la incertidumbre que se puede generar en 
la estimación de k y A a causa de la medición de las variables que 
contiene el modelo logarítmico, como la velocidad media u, la 
coordenada vertical y, la viscosidad cinemática v y la velocidad de corte 
u,. Parte de este análisis fue desarrollado por Segalini et a/. (2013). 


Conclusión 
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El objetivo de este documento fue estimar los valores x y A del modelo 
logarítmico a partir de datos experimentales medidos en un canal por 
los métodos de mínimos cuadrados, función de indicador, dispersión 
xk —A y topología. 


El método propuesto de topología permite estimar de forma simultánea 
k y A, ubicando el valor óptimo de E,,;, en un plano de estimación de 
error. De acuerdo con los valores obtenidos de x y A, y con la evaluación 
de los errores Eja Y Eca, con el método de topología se logran mejores 
resultados, en comparación con los métodos de función de indicador y 
dispersión kx — A, y con tendencia similar comparando con el método de 
mínimos cuadrados utilizando un valor de norma p = 2. Para el caso de 
norma p = 1 se tiene una mejor estimación en el método de topología en 
comparación con mínimos cuadrados, utilizando como criterio E, . 


Además, se comprobó el ajuste de la ley logarítmica con los valores 
obtenidos de la prueba experimental en el rango 30 < y* <0.36*, 
utilizando los valores calculados de x y A con el método de topología. 
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